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We investigate gluon propagators and the effective mass of the gluon fields in the MA gauge with 
U(l)3xU(l)8 Landau gauge fixing in SU(3) lattice QCD. The Monte Carlo simulation is per- 
formed on 16^ at )3=5.7, 5.8 and 6.0 and 32^ at j3 =5.8 and 6.0 at the quenched level. To calculate 
the propagators, we adopt a method to extract gauge fields from link-variables analytically in the 
SU(3) case. The off-diagonal gluons behave as massive vector bosons with the approximate ef- 
fective mass Moff ~ 1.1 — 1 .2Ge V in the region of r = 0.3 — 0.8fm, and the propagation is limited 
within a short range. On the other hand, the diagonal gluons behave as light vector bosons with 
A^diag — 0.3GeV and the propagation of diagonal gluons remains even in a large range. In this 
way, infrared Abeiian dominance is shown in terms of short-range propagation of off-diagonal 
gluons. Furthermore, we investigate the functional form of the off-diagonal gluon propagator 
The functional form is well described by the four-dimensional Euclidean Yukawa-type function 
exp(— moff'") / with m^^ =1.3—1 .4Ge V for r = 0. 1 — 0.8fm. This also indicates that the spectral 
function of off-diagonal gluons has the negative-value region. 
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1. Introduction 

For the quark-confinement mechanism, the dual-superconductor picture was suggested by 
Nambu, 't Hooft and Mandelstam In this picture, there occurs color-magnetic monopole con- 
densation, and then the color-electric flux between the quark and the antiquark is squeezed as a one- 
dimensional tube due to the dual Higgs mechanism. From the viewpoint of the dual-superconductor 
picture in QCD, however, there are two assumptions of Abelian dominance [§, ||] and monopole 
condensation. Here, Abelian dominance means that only the diagonal gluon component seems to 
be significant to confinement. 

The various lattice QCD Monte Carlo simulations support these assumptions when the maxi- 
mally abelian (MA) gauge fixing is peifomied |§ g § ^ |, |, jT^, [Tl], 0, |T|]. 



According to these studies, only the diagonal gluons play a dominant role for the infrared 
QCD physics, which is called "infrared Abelian dominance". Infrared Abelian dominance means 
that off-diagonal gluons are not significant to infrared QCD. Therefore, the essence of infrared 
Abelian dominance lie in the behavior of the off-diagonal gluon propagator. 

The gluon propagators in the MA gauge has been investigated in SU(2) lattice Monte Carlo 



simulations [ ]10| , |13| , |14| ]. To investigate the gluon propagators in the MA gauge, it is desired to ex- 
tract the gluons exactly from the link- variables, because the Unk- variable cannot be expanded even 
for a small lattice spacing due to large fluctuation of gluons. In SU(2) lattice case, the extraction 
is easy to be done without any approximation, because of the SU(2) property. With this extraction, 
the SU(2) lattice simulation suggests that the off-diagonal gluons do not propagate in the infrared 



region due to the effective mass M^s ~ 1.2GeV, while the diagonal gluon widely propagates [10|. 

In this paper, we propose a method to extract the gluons from the link-variable directly and 
generally in SU(3) lattice QCD, and to investigate the gluon propagators in the MA gauge. 



2. Formalism to extract gluon fields from link- variables 



In this section, we consider a useful and general method to extract the gauge fields analytically 



and exactly from the link- variables whether |agA^(;c)| ^ 1 is satisfied or not 1 15]. 
To this end, we first define the hermite matrix, 

A = j.[U-U^)= ^.{e'"^^ - e-'"^^) = sinagA. 



For simplicity, we have omitted the Lorentz index and space-time arguments. 

Arbitrary hemiite matrix A can be diagonalized by a unitary transformation as 



/A 



Arf = nAn^ 







A3 



where D gSU(3). We can obtain the eigenvalues A,- (/ = 1,2,3) by solving 

det(xl - A) =0. 



(2.1) 



(2.2) 



(2.3) 
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This is a cubic equation on x with real coefficients. The eigenvalues A,, i.e., the solutions of the 
equation are 



where 



and 



^o,± =zo,±Va^ + i3/3 + a, (2.4) 



a = ^TrA G R, jS = -(A22A33 + A33A11 + AnA22 - IA23P - IA31P - |Ai2|^) G R. (2.6) 

Here, = 0(a, jS) G R is analytically obtained and the derivation is given by Cardano's method. 
In this way. A; is obtained. 

The unitary matrix Q. can be also derived as follows. By solving Ae; = Ajg,, we obtain normal- 
ized eigenvectors e,- = '{xi,yi,Zi) {i = 1,2,3). Because of \e\ = 1, there is non-zero component, so 
that Zi is assumed to be nonzero without loss of generality and rescale it by l/zi, 

A ' {xi/zi yi/zi 1 ) = A; ' {xi/zi yi/zi 1 ) . (2.7) 

This is solved easily as 

From the normalization condition = 1, we obtain Q.^ = {ei,e2,es). 

When we diagonalize A with the unitary matrix the gluon fields A are diagonalized. Thus, 
the link-variables U = e'"^^ are also diagonalized with the unitary matrix Q. as 

Ud = QUQ^ = e'"^^^' = diag f e'^' , e'^^ , e'^^ ) , (2.9) 



where — tt < 0; < TT (/ = 1,2,3) is taken. 

In this way, we can derive gluon fields A from link-variables U analytically by diagonalizing 
them. 



(2.10) 



This formalism is quite general, because the derivation is correct in any gauge and even without 
any gauge fixing. 

3. SU(3) lattice QCD results for gluon propagators in the MA gauge 




Using the SU(3) lattice QCD, we calculate the gluon propagators [16] in the MA gauge with 



the U(l)3 xU(l)8 Landau gauge fixing. In the MA gauge, to investigate the gluon propagators, we 
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use the above-mentioned method. The Monte Carlo simulation is performed with the standard pla- 
quette action on the 16^ lattice with /3 =5.7, 5.8 and 6.0 and 32^ with /3 = 5.8 and 6.0 at the quenched 
level. All measurements are done every 500 sweeps after a thermalization of 10,000 sweeps using 
the pseudo heat-bath algorithm. We prepare 50 configurations with 16^ and 20 configurations with 
32"^ at each j8. The statistical error is estimated with the jackknife method. 

Here, we study the Euclidean scalar combination of the diagonal (Abelian) and off-diagonal 
gluon propagators as 

Gf;!\r) = ^ I Wx)Al{y)) , Gf^{r) ^ ^ £ {Al{x)A^,{y)) . (3.1) 

The scalar combination of the propagator is expressed as the function of the four-dimensional 
Euclidean distance r = {x^ —y^)^- When we consider the renormalization, these propagators 
are multiphed by an r-independent constant, according to a constant renormalization factor of the 
renormalized gluon fields. 

We show in Fig.|l| the lattice QCD result for the diagonal gluon propagator G^^^{r) and the off- 
diagonal gluon propagator G°^^{r) in the MA gauge with the U(l)3xU(l)8 Landau gauge fixing. 
In the MA gauge, G^^^{r) and G°^^{r) manifestly differ. The diagonal-gluon propagator G^^^{r) 
takes a large value even at the long distance. In fact, the diagonal gluons in the MA gauge 

propagate over the long distance. In contrast, the off-diagonal gluon propagator G^Jj^ (r) rapidly 
decreases and is negligible for r ^ 0.4fm in comparison with G^^^\r). Then, the off-diagonal 
gluons A° {a ^ 3,8) seem to propagate only within the short range as r ;^ 0.4fm. Thus, "infrai^ed 
abelian dominance" is found in the MA gauge. 



a 1.5 





Figure 1: The SU(3) lattice QCD results of the gluon propagators G^j^\r) and G°^(r) as the function of 
r = (x^ ~yiJ.y in the MA gauge with the U(l)3xU(l)8 Landau gauge fixing in the physical unit. The 
Monte Carlo simulation is performed on the 16^ lattice with j3 = 5.7, 5.8 and 6.0. The diagonal-gluon 
propagator G^Jf'(/") takes a large value even at the long distance. On the other hand, the off-diagonal gluon 
propagator G?E (r) rapidly decreases. 



4. Estimation of diagonal and off-diagonal gluon mass in the MA gauge 



Next, we investigate the effective mass of diagonal and off -diagonal gluons [16|. We start from 
the Lagrangian of the free massive vector field with the mass M / in the Euclidean metric. In 
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the infrared region with large Mr, the propagator G^^{r,M) reduces to 

3 M 3\/M e^'^'' 

G,,ir,M) = {A,ix)A,iy)) = 7^-^i(Mr) ^ -f-y^, (4-1) 

^Tt r 2(271)2 r2 




r[fm] 



Figure 2: The logarithmic plot of r^l'^G°^^{r) and r^^^G^^^{r) as the function of the Euclidean distance r 
in the MA gauge with the U(l)3xU(l)8 Landau gauge fixing, in the SU(3) lattice QCD with 16"* at )3 = 5.7, 
5.8 and 6.0. The sohd line denotes the logarithmic plot of r^l^G^^i (r) ~ r'/^A^i (Mr) in the Proca formalism. 




Figure 3: The logarithmic plot of r^^^G°Jj^{r) and r^^^G^^^{r) as the function of the Euclidean distance r 
in the MA gauge with the U(l)3xU(l)8 Landau gauge fixing, in the SU(3) lattice QCD with 32^* at j3 = 5.8 
and 6.0. 



In Fig.^ we show the logarithmic plot of r^^^G"Jj^{r) and r^^^G^^^{r) as the function of the 
Euclidean distance r in the MA gauge with the U(l)3 xU(l)8 Landau gauge fixing. From the linear 
slope on r^^^G°Jj^{r) in the range of r = 0.3 — 0.8 fm, the effective off-diagonal gluon mass Moff is 
estimated. Note that the gluon-field renormalization does not affect the gluon mass estimate, since 
it gives only an overall constant factor for the propagator. We summarize in Table 1 the effective 
off-diagonal gluon mass Moff obtained from the slope analysis with 16^ at jS =5.7, 5.8 and 6.0. 
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Therefore, the off-diagonal gluons seem to have a large mass Moff ~ 1.1 — 1.2 GeV. This result 
approximately coincides with SU(2) lattice calculation [|lO|]. 

On the other hand, for the diagonal gluons, their propagator seems to have some dependence 
on j8 and lattice volume. Therefore, we estimate the effective diagonal gluon mass M^i^g with larger 
lattice size, 32^*. In Fig.^, we present the logarithmic plot of r^^^G°^^{r) and r^^^G'lf^^{r) with 32"^ 
at j8 =5.8 and 6.0. From the linear slope on r^^^G'j^^^{r) in the range of r = 0.3 — 0.8 fm, the 
effective diagonal gluon mass is estimated as Mdiag — 0.3GeV at each j8. 

Table 1: Summary table of conditions and results in SU(3) lattice QCD. In the MA gauge, the off-diagonal 
gluons seem to have a large effective mass Moff ~ 1.1 — 1.2 GeV and the functional form in the range of 
r = 0.1 — 0.8 fm is well described with the four-dimensional Euclidean Yukawa function ^ exp(— OToffr)/r 
with moff ~ 1 .3 - 1 .4 GeV. 



lattice size 


P 


fl[fm] 


Moff[GeV] 


moff [GeV] 




5.7 


0.186 


1.2 


1.3 


164 


5.8 


0.152 


1.1 


1.3 




6.0 


0.104 


1.1 


1.4 



Finally in this section, we discuss the relation between infrared abelian dominance and the 
off-diagonal gluon mass. Due to the large effective mass Moff, the off-diagonal gluon propagation 
is restricted within about M^^ ~ 0.2fm in the MA gauge. Therefore, at the infrared scale as r » 
0.2fm, the off-diagonal gluons {a ^ 3, 8) cannot mediate the long-range force like the massive 
weak bosons in the Weinberg-Salam model, and only the diagonal gluons A^, can mediate the 
long-range interaction in the MA gauge. In fact, in the MA gauge, the off-diagonal gluons are 
expected to be inactive due to the large mass Moff in the infrared region in comparison with the 
diagonal gluons. Then, infrared abelian dominance holds for r S> M^^^jf . 

5. Analysis of the functional form of off-diagonal gluon propagator in the MA gauge 

In this section, we investigate the functional form of the off-diagonal gluon propagator G^^Jj (r) 



in the MA gauge in SU(3) lattice QCD [16]. In the previous section, we compare the off-diagonal 
gluon propagator with the massive vector boson propagator and estimate the gluon mass. In fact, 
the gluon propagator would not be described by a simple massive propagator in the whole region 
of r = 0.1-0.8fm. 

There is the similar situation in the Landau gauge [|l7]|. The functional form of the gluon 
propagator cannot be described by exp(— Mr)/r^/^ with an effective mass M in the whole region 
of r = 0.1 — 1.0 fm. The appropriate form is the four-dimensional Euclidean Yukawa-type function 
exp(— mr)/r with a mass parameter m. 

In the same way, in the MA gauge, we also compare the gluon propagator with the four- 
dimensional Euclidean Yukawa function. In Fig.^ we show the logarithmic plot of ''G^^ (r) and 
rG^I^'^'(r) as the function of the distance r in the MA gauge with the U(l)3xU(l)8 Landau gauge 
fixing. Note that the logarithmic plot of rG°^^{r) is almost linear- in the whole region of r = 0.1 — 
0.8 fm, and therefore the off-diagonal gluon propagator is well expressed by the four-dimensional 
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Euclidean Yukawa function Ae '""'f' /r, with a mass parameter nton and a dimensionless constant 
A. The best-fit mass parameter nion is given in Table 1 at each j8 = 5.7, 5.8 and 6.0. 



V A 



0.4 0.5 
r [fm] 



P=6.0 
P=6.0 
P=5.8 
P=5.8 
P=5.7 
P=5.7 



Figure 4: The logarithmic plot of /"GJ^^ (r) and rG^^^ (r) as the function of the Euclidean distance r in the 
MA gauge with the U(l)3xU(l)8 Landau gauge fixing, using the SU(3) lattice QCD with 16'^ at j3=5.7, 5.8 
and 6.0. For rG^^^ (r), the approximate hnear correlaton is found. 



We comment on the four-dimensional Euclidean Yukawa-type propagator [17|. If the func 



tional form of the off-diagonal gluon is well described by the four-dimensional Yukawa function, 

-off 



we analytically calculate the off-diagonal zero-spatial-momentum propagator, D^^{t) = J d^xG°Ji{r), 



and obtain the spectral function p°^{(o) hy the inverse Laplace transformation in the MA gauge: 



p°"{co) = -— TTJn^i^ - "^off + 7- uTI^ « - '"off ■ (5.1) 

{(O^-m^^y/^ (w-moff)^/^ 

6. Summary and Concluding Remarks 

We have performed the study of the gluon propagators in the MA gauge with the U(l)3 xU(l)8 
Landau gauge fixing in the SU(3) quenched lattice QCD. To investigate the gluon propagators in 
the MA gauge, we have considered to derive the gluon fields analytically from the SU(3) link- 
variables. 

With this method, we have calculated the Euclidean scalar combination G^^(r) of the diag- 
onal and the off-diagonal gluon propagators, and have considered the origin of infrared Abelian 
dominance. The Monte Carlo simulation is performed on the 16'^ lattice at j8=5.7, 5.8 and 6.0 and 
on the 32^ at j8=5.8 and 6.0 at the quenched level. We have found that the off-diagonal gluons be- 
have as massive vector bosons with the effective mass Mos ~ 1.1 — 1.2 GeV for r = 0.3 — 0.8 fm. 
The effective gluon mass has been estimated from the linear fit analysis of the logarithmic plot of 
r^/^G^^ (r). Due to the large value, the finite-size effect for the off-diagonal gluon mass is expected 
to be ignored. The large gluon mass shows that the off-diagonal gluons cannot mediate the interac- 
tion over the large distance as r » M^^ , and such an infrared inactivity of the off-diagonal gluons 
would lead infrared Abelian dominance in the MA gauge. 
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On the other hand, from the behavior of the diagonal gluon propagator G^^^{r) and r'^^G^^^{r), 
the diagonal gluons seem to behave as light vector bosons with Mjiag ^0.3 GeV for r = 0.3 — 0.8 fm 



[16], considering also the larger- volume analysis with 32^ at /3 =5.8 and 6.0. 

Finally, we have also investigated the functional form of the off-diagonal gluon propagator 
G^^(r) in the MA gauge. We show that G^^(r) is well described by the four-dimensional Eu- 
cUdean Yukawa-type form with the mass parameter /Woff ^ 1.3 — 1.4 GeV in the whole region of 



r = 0. 1 — 0.8 fm. This indicates that the spectral function p°^(co) of the off-diagonal gluons in the 



MA gauge has the negative- value region [[161], as in the Landau gauge [|1^, |T8|, [T9| ]. 
Acknowledgements 

This work is supported in part by the Grant for Scientific Research [(C) No. 23540306, Priority 
Areas "New Hadrons" (£01:21105006)], Grant-in-Aid for JSPS Fellows (No.23-752, 24-1458) 
from the Ministry of Education, Culture, Science and Technology (MEXT) of Japan, and the Global 
COE Program at Kyoto University. The lattice QCD calculations are done on NEC SX-8R at Osaka 
University. 

References 

[1] Y. Nambu, Phys. Rev. DIO (1974) 4262; G. 't Hooft, in High Energy Physics (1975); 
S. Mandelstam, Phys. Kept. 23 (1976) 245 . 

[2] G. 't Hooft, Nucl. Phys. B190 (1981) 455. 

[3] Z.F. Ezawa and A. Iwazaki, Phys. Rev. D25 (1982) 26U\Phys. Rev D26 (1982) 631. 

[4] A. S. Kronfeld, G. Schierholz and U.-J. Wiese, Nucl. Phys. B293 (1987) 461; 

A. S. Kronfeld, M. L. Laursen, G. Schierholz and U.-J. Wiese, Phys. Lett. B198 (1987) 516. 

[5] T. Suzuki and I. Yotsuyanagi, Phys. Rev D42 (1990) 4257. 

[6] F Brandstaeter, U.-J. Wiese and G. Schierholz, Phys. Lett, mil (1991) 319. 

[7] J. D. Stack, S. D. Neiman and R. J. Wensley, Phys. Rev D50 (1994) 3399. 

[8] O. Miyamura, Phys. Lett. B353 (1995) 91. 

[9] R. M. Woloshyn, Phys. Rev D51 (1995) 641 1. 

[10] K. Amemiya and H. Suganuma, Phys. Rev D60 (1999) 1 14509. 

[11] H. Ichie and H. Suganuma, Nucl. Phys. B548 (1999) 365; Nucl. Phys. B574 (2000) 70. 

[12] H. Suganuma et al., Nucl. Phys. Proc. Suppl. 106 (2002) 679; A. Shibata et al., PoS (Lattice2007) 331. 

[13] V. G. Bornyakov, M. N. Chernodub, F. V. Gubarev, S. M. Morozov and M. I. Polikarpov, 
Phys. Lett. B559 (2003) 214. 

[14] A. Cucchieri, F Karsch and R Petreczky, Phys. Lett. B497 (2001) 80. 

[15] S. Furui and H. Nakajima, Phys. Rev D69 (2004) 074505; Nucl. Phys. Proc. Suppl. 73 (1999) 865. 
[16] S. Gongyo, T. Iritani and H. Suganuma, Phys. Rev D86 (2012) 094018; PoS (Lattice2012) 212. 
[17] T. Iritani, H. Suganuma and H. lida, Phys. Rev D80 (2009) 1 14505. 
[18] J. E. Mandula and M. Ogilvie, Phys. Lett. B185 (1987) 127. 

[19] R O. Bowman et al., Phys. Rev D70 (2004) 034509; Phys. Rev D76 (2007) 094505. 



8 



